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Abstract 

o 



We study the linear perturbation of the recently proposed model of inflation where a uniform 
gauge-kinetic coupling of the inflaton to multiple vector fields breaks the cosmic no-hair conjecture 
while maintaining the isotropy. We derive the general quadratic action for the perturbation and 

u 

calculate the power spectra of scalar and tensor modes at the end of inflation by in-in formalism. 
It is shown that the model predicts slightly red spectra and the tensor-to-scalar ratio tends to be 

o: 

suppressed. The comparison with the data from WMAP 7-year does not impose strong constraints 
on the parameters and both weak- and strong- gauge-field regimes are consistent with the current 



observations. 



* K.Yamamoto@damtp.cam.ac.uk 



I. INTRODUCTION 



Over the past decade, the prospect of unveiling minute details of density fluctuations 
through observations of the Cosmic Microwave Background (CMB) and large scale struc- 
ture has driven a growing interest in inflationary scenarios beyond the single-field slow-roll 
model. Despite its simplicity and the remarkable success in matching the observed density 
power spectrum, it is hard to believe that a single scalar field is entirely responsible for the 
dynamics of the early universe. Along with multi-scalar models such as hybrid inflation 1], 
attempts have been made to incorporate vector fields 2j-|4|, which have proven to be largely 
unsuccessful due to various instabilities in realizing accelerated expansion p-l9]. Recently, a 
viable model of vector inflation has been proposed by Soda and his collaborators 10] where 
the gauge-kinetic coupling with inflaton, originating from supergravity, maintains the am- 
plitude of the vector and results in non-trivial signatures in the primordial fluctuations, e.g. 
statistical anisotropy and scalar-tensor correlations 11] . This "inflation with vector-hair" 
has been subsequently scrutinized 12Ml71] and its stability has been widely established \18h 
201 ] . Moreover, it was found that the inclusion of additional vector fields tends to reduce 
anisotropy in the background through a generic dynamical mechanism 21] . In the case 
of uniform gauge-kinetic coupling for three or more vectors, the trajectories converge to a 
universal isotropic attractor with non-vanishing vector energy density, whose amplitude is 
determined by the strength of the coupling. Very similar dynamical behaviors have been 
observed in the models employing non-Abelian gauge fields 22H25] . which indicates that 
the tendency towards isotropy is generic in the inflation with multiple vectorial degrees of 
freedom. This special case is indistinguishable from the single-field slow-roll model at the 
level of classical background dynamics as long as it provides a sufficient e-folding number. 
As such, it is necessary to investigate its linear perturbation in order to impose observational 
constraints and decide its viability, which is the aim of this article. 

Incidentally, these gauge-kinetic models can be viewed as the classical counterpart for the 
particle production effects in preheating scenarios. In this context, couplings between the 
inflaton and gauge fields have often been discussed in the attempts of generating primordial 
magnetic fields |26l . 127] . gravitational waves [28], and non-Gaussianity [29] . The crucial 



conceptual leap of the inflationary models considered in this article is the breach of the 
cosmic-no-hair conjecture whereby typical preheating scenarios assume vanishing classical 



background values of the fields excited by the inflaton, which forces the analysis to go 
beyond linear order and makes it complicated. The presence of vector-hair in inflation 
illustrates that the non-vanishing background energy density of (potentially anisotropic) 
auxiliary fields does not necessarily mean a break down of the inflationary regime and the 
accelerated expansion may continue without wiping out classical "hair," being the attractor 
solution in the phase space at the same time. Therefore, it is interesting to ask whether the 
existence of the background gauge fields can produce any of the features mentioned above 
within the well-established linear perturbation around the isotropic Friedmann-Lemaitre- 
Robertson- Walker space-time and quantum field theory in the quasi-de-Sitter background, 
and how large is the amplitude if any. 

In the present article, we calculate the spectra of the scalar curvature in uniform density 
slicing and the gravitational waves generated through the interaction of the perturbative 
variables with the background triad of the gauge fields. The result obtained is very much 
analogous to what was found for the anisotropic cases [ill, |3C], |3l( . We find that the scalar 
power spectrum acquires a double logarithmic scale dependence with the magnitude of the 
correction being given by the fractional energy density of the background gauge fields with 
respect to the scalar kinetic energy. The effect is not slow-roll suppressed due to the steepness 
of the gauge-kinetic function that is needed to maintain the background gauge fields in the 
accelerated expansion. While the tensor mode receives similar corrections, the magnitude is 
smaller than that for the scalar mode by a factor of slow-roll parameter, which was also seen 
in the anisotropic cases. However, since the effect is completely isotropic, the restrictions 
coming from the observational data are much weaker than the anisotropic single-gauge-field 
models. Although the correction terms involve e-folding numbers, which resulted in strong 
constraints on the energy density of the gauge field in the anisotropic models, the large face- 
vale of the e-folding numbers can only affect the overall amplitude of the power spectra in 
this isotropic setup. Since the amplitudes of the quantum fluctuations are always normalized 
by the energy scale of inflation, the potentially large corrections to these amplitudes can be 
absorbed into this normalization, whence do not immediately pose a problem. 

The paper is organized as follows. In section II, we present the background equations 
and introduce the slow-roll parameters to characterize the inflationary dynamics. Section 
III explains the structure of the second order action for the gauge fields, which lead to an 
extended notion of scalar-vector-tensor decomposition. In section IV, the power spectrum 



of the scalar mode is computed by employing in-in formalism on the de-Sitter space-time. 
In section V, a similar calculation is performed for the tensor mode and observables such as 
spectral tilt and tensor-to-scalar ratio are obtained. Section VI summarizes the results and 
discusses the outlook. 

II. THE ISOTROPIC BACKGROUND DYNAMICS OF THE INFLATION WITH 
GAUGE-KINETIC COUPLING 

We consider the universe described by the following action; 

S = fdx^^g ^&R - \d^d^ - V{ip) - ^ ( ;l F^A . (1) 

The scalar field (p acts as the inflaton and interacts with the gauge fields through the gauge- 
kinetic function f((p). Although we assume for simplicity that Ffi?\ (m = 1,2,3) are three 

copies of Abelian gauge field 

F M =(cM M)^ (2 ) 

there is strong evidence for the fact that the dynamics of an 577(2) gauge field with a 



gauge- kinetic coupling in an inflationary regime is well described by the action ([T]) 15] . The 
generalization for m > 4 is also straightforward. We adopt the ADM formalism and follow 
the notation of the ref. |32| to write the metric as 

I -iV 2 + N k N k N 3 \ 

g»v = Ar ( 3 ) 

where 

g ik g kj = S^, N* = g^N y (4) 

It is well known that in terms of the normalized extrinsic curvature 

Ea = -\{9a-2N m ) (5) 

and the intrinsic scalar curvature of the constant time slice 

3 R = (gi jM + j»r™r»J (s V - g ij g kl ) , (6) 

the gravitational part of the Lagrangian is written as 

^9 = : §faE*-E?) + ±NVg*R. (7) 



We introduce the electric fields by 

E M = F M (g) 

for the 1 + 3 decomposition and write the gauge-field Lagrangian in the form 

Cm = ^f 2 9 lk (E^ + F^N*) (e™ + F^N 1 ) - ^A^VV^^]- (9) 
As usual, the scalar Lagrangian is given by 

C v = & Q^ 2 - W,iN l + fairfA - N^ Q<7>,^ + %)j . (10) 

For the background, we use the ansatz 
N=Af(t), N t = 0, 9ij = a(t)%, <p = <p(t), A { ™ ] = 0, A { ™ ] = A{t)5 m % . (11) 
The resulting equations of motion read 

SM^-^-ATV-gi^O, (13) 

$ + 3~0 ~ §- r + N*V„ - ^%i 2 = 0, (14) 

a TV ar 

^A = const. = cM p i. (15) 

The first integral (fT5l) will be used to eliminate A. In order to figure out the conditions 
for inflation, let t be the proper time coordinate by choosing M = 1 and define the Hubble 
expansion rate H = a/a. For the spatial slice to undergo accelerated expansion and for it 
to last more than one Hubble time, it is required that the parameters 

= _H_ = 2 . c 2 , 1fi * 

tR H 2 2M 2 l H 2 a A f 2 H 2 ' l ' 

''«-7fe- 2 (- + 2^i|)- 2 ( 1 -2^)( 2+ 7rf) < 17 > 

be much less than unity These two parameters characterize the evolution of the space-time: 
e# <C 1 guarantees accelerated expansion a/ a > and constancy of H over a few Hubble 
times and ^ < 1 ensures that this regime lasts for at least tens of e-foldings. Another 



important element is the balance between the scalar kinetic energy and the amplitude of the 
gauge fields. It proves to be convenient to use the following set of parameters; 

** = l^p < ^ ( 18 ) 



*=^ = H e ' + ^;- (l9) 

e v represents the kinetic energy of inflaton, which is already much smaller than unity by the 
conditions above. Its significance resides in the fact that it controls the balance between the 
kinetic energy of the inflaton and the gauge fields through 



f 2 A 2 c 2 M 2 , 



a 2 M^H 2 a 4 f 2 H 2 



e H -e v < e H . (20) 



In connection to this balance, it later proves to be useful to define another parameter 



' ch — £ v 



1=J^-^, (21) 

e f 

measuring the ratio of kinetic energy between the background gauge fields and the inflaton. 

Note that there is no a priori constraint on X as long as e v < e# <C 1 is satisfied. Although it 

would not be necessary to assume ^ < 1, a large rj^ means rapidly varying e v and in turn a 

rapid exchange of the energy between the gauge fields and the inflaton, which would require 

a dedicated condition to be met by V and /. In order to avoid unnecessary complications, 

we do assume the smallness of r] v , which leads to the control over the gradients of V and / 

through 

V v _ 1 / .. ., ., 3 1 



M pl H 2 JTe 



[ -6e H + 3e 2 H - -e H i] H - e H e^ + -e^ J (22) 



M piU _ l ( o , l enVH-e^ . 



and 

to 

f y/2e^ V 2 e H ~ ^ 

The slower (p rolls down the potential, the greater is the amplitude of the gauge fields and 
the steeper is the slope of the gauge-kinetic function. 

III. THE SECOND ORDER ACTION AND SCALAR- VECTOR-TENSOR DE- 
COMPOSITION 

Now let us perturb the background and write down the second order action. Our con- 
vention for the perturbative variables are given as follows: 
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Metric: 



Inflaton: 



Gauge fields: 



N = M(l + <f>), Ni=MaPi, g ij = a 2 (6 ij + 2 lij ). (24) 



ip = (p + TV. (25) 



4 m) = a m , A M = m + x™. (26) 



Note that the definition of <fi is unconventional. The derivation for the gravity and scalar 



actions can be seen in any standard literature (e.g. ref. 33] ). To derive the gauge- field 



Lagrangian, it is convenient to define the perturbed electric fields 

XI = n - v™. (27) 

Substituting 

E (m) = ^ra + ^ f H = _ 2 ^ ^ (2g) 

and 

E (m) + F (rn) NJ = ^ + X™ - 2— X \^ j + higher order terms, (29) 

into the action (|9]), one obtains 

~n^U - 2/^xWi + ^V " w ^ 7 + 6// -^ (30) 

2Af V2 W V + /V 7+ 4 AT 

where 7 = 7^. Combining this with the gravity and scalar Lagrangians, using the back- 
ground equations and performing integration by parts, we derive the general quadratic 
Lagrangian as 

(2 a 3 M* NaMl 

O > = -^- [liffij - 7 ) + — - — (27ij,i7ifc,fc - lij,klij,k + ^llij.H ~ 11, a) 

+ 2M A ^^ - 7 ) + 2M n ~ TW - 2 T a Kw A?~ A J ff 

a/ 2 . . Aff 2 k k gf ■ 2a//,,, . 

27V J J ' ~ "^T* [M ' ]X M ~ AT 7ij j ~A/^ m 



+^H^ OWfci) - Pli) + « 2m piA ( lijj - 7,i + 2^0, 



pi / ~ „o \ o , ,0 ^ / . . .. ^ 



a 



-Afa 3 V(p 2 + 7 



a 3 (p 

■7T 



-a 2 (3i ( 0ir,i + 2±-A x \ kA 



P 

a 



Afa?V 



:<P 



+ M^ 



A/" 
2a 2 a 



^ 



7T 



— 4)f m 



+<t> 



af 



i2 



A^ 



a 3 cp . 

-TV 



Na*V M + 






7 + Ma ('jijjj - 7^) 
a/ 2 



(31) 



7T 



4 Y m 

A" m 



AT ' AT \J * ' AT 

Although the upper indices originated from the label attached to each copy of gauge field, 
some of them are now contracted with spatial indices. This happened since the background 
respects the spatial 0(3) symmetry and is also invariant under linear mixing of the three 
vector fields by 0(3) "rotation," namely A \ oc 5™, which resulted in transferring spatial 
indices to the ones denoting species. Consequently, x % j ( an d X l j) looks as if it were a 3 x 3 
spatial tensor even though it represents a triplet of vectors. This suggests that one can 
formally extend the scalar-vector-tensor decomposition of the perturbed quantities as was 
first carried out in the ref. [22J. We use the following symbols; 






B a — Si 



1 



-ip5ij + E ti j + F(jj) + -hij, 



X j — a $ij + 0,ij + e ijk T ,k + K (i,j) + tijk^k + <^ij i 
a 1 = (Aj + Vi. 



(32) 
(33) 
(34) 
(35) 



As usual, Si,F t ,Ki, Aj and z/j are divergence-free vectors and h^ and Wy are traceless trans- 
verse tensors and e^ is the Levi-Civita symbol in three dimensions. Now all the indices 
are downstairs, which signals we are going to ignore the distinction between spatial and 
component indices. It can be easily seen that the Lagrangian indeed splits into three pieces 
each of which contains only scalars, vectors and tensors respectively, up to a surface term, 
and one can deal with each mode separately as far as linear perturbations are concerned. 

Not all of the quantities are dynamical because of the gauge freedom of general relativity 
and the U(l) gauge symmetry. Let us first consider an infinitesimal space-time diffeomor- 
phism 

t-tt + ri, x l -> x i + + &■ (36) 

It induces the following transformations: 



fj, -)• fji + A£, Vi-* Ui + A& 
a-ta-Ar), 9 -> 9 + A£ 



(37) 
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Note that r as well as u;^- are gauge invariant. One is also allowed to perform the infinitesimal 
U(l) transformations 

where p^ are arbitrary scalar functions. It is again convenient to decompose them as 

P {m) = P,m + Pm (39) 

and the transformation laws become 

P^-P + P, Ui^Ui + Pi, 

9 ^6 + p, Ki-> Ki + pi, (40) 

Aj — > Aj + —tijkPj,k- 

It is instructive to count the numbers of degrees of freedom. There are four scalars from 
the metric, four from the gauge fields and another for the inflaton. Two from the metric 
(<f), B) and one from the gauge fields (p) are non-dynamical. We have three arbitrary func- 
tions for gauge transformations, which results in three dynamical scalar degrees of freedom. 
For vector perturbations, there are two from the metric and three from the gauge-fields, 
among which we have one non-dynamical for each sector (Si and v,j). Two more are redun- 
dant because of the gauge freedom and we are left with only one dynamical vector mode. 
Finally, there are two gauge-invariant dynamical tensor degrees of freedom. 

IV. SCALAR MODES AND THE CURVATURE POWER SPECTRUM 

In this section, we consider the scalar perturbation; 0, if), B, E, n, a, 6, r and p and com- 
pute the curvature power spectrum. We fix the gauge by setting 

iP = E = 9 + p = 0, (41) 

which means the constant-time hypersurfaces are flat and n, a and r are the three dynamical 
variables. The convenient choice of time-coordinate is the conformal time r] (M = a) and 
differentiation with respect to it is denoted by primes. We set M v \ = 1 in this section and 
the next to save the space. 
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A. The Lagrangian and the Curvature of the Uniform-density Hypersurfaces 



The second order Lagrangian is written as 

4 2) = y (** " ™r,0 - 5 ( A w " 3 (//,„)„ A' 2 ) vr 2 + 6//,^W 

+/ 2 (Jc/ 2 - a,ia,i + r'^ - TuTjA - a 2 B U%,u - <f\u - 2 ^ A ' a A ( 42 ) 
-aV0 2 - [a^V + (aV^ + Zff^A' 2 ) tt + 3/ 2 AV] . 

Varying B gives 

6 = —Ti + -!—=A!a (43) 



a 1 



where we defined 

H = -. (44) 

a 

Its substitution into the action eliminates the non-dynamical variables <fi and B and leaves 
three scalar degrees of freedom. Using the canonically normalized variables 

an, 

a = VSfa, (45) 

f = V2fr, 

and performing several integrations by parts, it yields 

4 2) = \ 0' 2 - (Vvr) 2 - m 2 7r 2 ) + \ (ol 2 - | (Vd) 2 - m 2 d 2 ) 

1 / f" 

I /V7^V2 /w2^\2 i J /V7A\2 

where 



+ 2( v (Vr) /2 -(V 2 f) 2 + y(Vf 



,2 



^„ 



-(2 + 3e v - e H - e^e* + r) (p e ip )'H 2 



(V w + v^V; + 6^) - 3(6* - e„) (^ + jf - V^y) « 2 , 



f" 2 

,2 J , /, , \ /o , \ 1/2 z /, , \ „2i 



< = -y + (en - e^) (3 - e H ) U 2 - - (e H - e v ) a 2 V, (47) 

2 ^ (Cg 3 ~^ (?/ 2 - aV) - ^3(6* - e v ) (2^ - ^1) ^ 
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We notice that f is decoupled from the other degrees of freedom since it is a pseudo-scalar 
representing the magnetic fields, which are zero in the background and hence gauge invariant. 
As its evolution is trivial and only affects the curvature perturbation at a non-linear order, 
we ignore this contribution. It should be noted, however, that it may well play an important 
role in the three-point or higher order correlation functions. On the other hand, a makes a 
physical contribution to the scalar curvature at linear order directly (see the equation (1491) ) 
and through the interaction with n, even though its origin is the vectorial gauge fields. It 
is not surprising since this d-mode represents the modulation of the background value of 
the gauge triad which already affects the background expansion, a is a superposition of 
three mutually orthogonal vector modes, which also explains why its propagation velocity 
is reduced from the speed of light to c 2 = 2/3. 

Our aim is to compute the curvature perturbation in the uniform density gauge, which 
is given by 

-C = ^ + H% (48) 



in the linear order according to Malik and Wands 34j. This quantity is not conserved during 
the inflationary period in the present model, even on the super-horizon scales because of the 
multi-field interaction and we will calculate its spectrum at the end of inflation quantum 
mechanically by tree-level perturbative expansion of the in-in formalism. Although it might 
be modified even after inflation by the entropy perturbations potentially generated by the 
gauge fields, we will not discuss it as the details depend on reheating. The definition of the 
total energy perturbation 5p here is that of Kodama and Sasaki |35j . At linear order, it is 
just the 0-0 component of the energy-momentum tensor and we have 

5p = -8T° (49) 



GT V 7 !"' + 



* a + 44(4_£, +vl 



2Ha 2 a 4 f 2 \f 1U 
3 x /2^m' + a' 3 



OC * C //o OC 



2^(1 - eJW 2 + 3(£« -e v )[ ! -f-<ff)H 2 + a 2 V v 



7T 



+a- 3 y /3(e H -e ip )H&' + er 3 ^-^ ("y« + (Se H - t v )H 2 \ a. 

So far, we haven't used the smallness of the four parameters, €^,£^,77^, and 77^. Keeping 
the leading order terms in these parameters and assuming 

i) H = 0(1)H VH , ^ = 0(l)Hr} v , (50) 
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which is reasonable, the curvature perturbation is expressed as 

(51) 



c _ ' l e * 



A . 4 + 6X 2 . yz^( „, 2 A 

it H it + -t / —X a a 

V V2 V V 



Zl-iaen y 2 
whose evolution is governed by the Lagrangian 

£ = V - i(W) 2 + ^-^-vr 2 + -a' 2 - Uva) 1 + -n 2 + ^-^rn' ^f"n. (52) 



i^/2 1/«-n2 ! + 6X 2 l A/2 l /VT ^ 2 l.o 2 V / 6X_, 4v/6X 

^ _ (V7r) 2 + 2 + ,2 _ Va 2 + 2 + _V_£ a / _ _^L_ 

2 2 T] 2 3 r] z r] r] z 



The steepness of the gauge-kinetic function compensates the suppressions coming from the 
small background density of the inflaton and the gauge fields as well as the flatness of the 
inflaton potential, hence results in the apperance of the parameter X defined in ( 12 ip , which 
is potentially of order unity. 

B. Calculation of the Power Spectrum 

In order to calculate the two-point function, we treat the terms involving X as small 
perturbations. Then, n and a are free particles in the de-Sitter background at the leading 
order. We impose the Bunch-Davies vacuum condition and write them as 

j^y 3 (u k (vWe^ + vTMale-^) , (53) 

a( V , x) = J -0^ (u k (v)b^* + u* k (v)bie- ik ~) . (54) 

In the above expressions and all the following, the field symbols denote the corresponding 
quantities in the interaction picture. The vacuum mode functions are given by 

- ( "» =M i -Q e " b ' (55) 



- ( ") = vm I 1 " ^,) e ~" M (56) 

with c 2 = 2/3 being the propagation speed of a, and the creation and annihilation operators 
satisfy 

K> 41 = t 6 P' 6 i] = ( 27r )^(P - <$> t a P' 6 J = °> etc - ( 57 ) 

The interaction Hamiltonian is given by 

HM = J d 3 x [~^ 2 + ~^-^« + -tf-* a ) ■ (58) 
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Note that we performed integration by parts for simplifying the calculations below. In 
interpreting the expression, normal ordering is understood in order to eliminate vacuum 
bubbles. In terms of the de-Sitter vacuum state defined by 

a k |0) = 6k |0) = 0, (59) 

the in-state is written as 

| in) = Texp ( -i I dfjH^f])) |0) (60) 

V J-oo(l-ie) J 

where the symbol Texp denotes the time-ordered exponential. While the lower limit of the 
time integral should physically be the beginning of inflation, it does not affect the following 
calculations as long as it is far enough in the past. The quadratic variation of the curvature 
perturbation during inflation is defined by 



(27r) 3 ,5(k + p)(in|4 2 (r/)|in) = J d 3 x e ik - x (in|C(r ? ,x)C(r ? , 0)|in) 

= /^(in|(C k + C k )(Cp + Ci p )|m) (61) 

where we used 

/rl 3 k 
^ {(Me~' lk - X + C k (r/) f e- k - x ) . (62) 

Comparing this with f loTj) . we need to compute 

(in| (g k (vW + 9k(v) a -v) [9p(v)a p + 9p(TJ)a-p) l in )> ( 63 ) 

(in| (j k ( V )b k + f k ( V )bl k ) (jMK + f P (v)bl p ) |in), (64) 

and 

(in| \g k (r])a k + gl{r])a ] _^j {j p {rj)b v + j*(v)bi p ) |in), (65) 

where g k (v) an d 3k(v) are linear combinations of Uk(r)), Uk(vY an( ^ their derivatives given by 

4 + 6X 2 
9k(v) = Mv) + Mv) (66) 

2 

Mv) = u'kiv) - ~Uk(v)- (67) 

rj 

Let us start from fl63l) . Expanding the time ordered exponential up to second order in X, 

we need to keep the following contributions: 

|^(r ? )| 2 (27r) 3 5(k + p) (68) 
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drnHjirn) [g k (vW + <7*fa)aL k ) g* p (v)al p \0) (69) 

-i(0\g k (v)a k (g p (v)a p + 9p(vW- P ) J d»7i-H>(»7i)|0> (70) 

/v rvi , v 

tfefr / dr^Hj^Hjim) [g k (v)a k + g* k (v)al u ) <£fa)at p |0) (71) 

-(^\gk(v) a kyg P (v) a p + 9p(v) a -p) drji I dr] 2 H I (r] 1 )H I (r] 2 )\0) (72) 

d^H^rn) (s fe a k + ^a f _ k J ((7 p a p + ^alpj / drj^H^)^). (73) 



The first term in fl58l) . which can be written in momentum space as 
6X 2 f ,, . 6X 2 /" d 3 A; 



r/ 2 



/o o O-L I CL rC / 2 ii2t *2 t t "\ / \ 

d X7T = - / — — - I M fc a k a_ k + 2|wjfc| a^a^ + u k a^a'^j , (74) 

contributes to the second and third lines only at the order X 2 . This modified mass effect 
reduces to 

m> + m = 24X 2 (2n) 3 5(k + p) f d Vl \j m (u k (vi) 2 gl(v) 2 ) , (75) 

where 3m means taking the imaginary part. The integration can be readily carried out by 
using the cosine- and sine-integrals that are defined by 

Ci(x)= / dt— = 7 + lnx + 0(x 2 ), (76) 

Si(x) = - I* dt— = - + 0(x). (77) 

J CO t * 

It is convenient to introduce x = kr] and the rational functions 



X 2 J X 2 



3(1 + 2X2) / 3(1 + 2X2) 
by which g* 2 is expressed as 



Q(*)= x [ 1+ X 2 - )■ ( 79 ) 



gl(v) 2 = ^(-P(x) + 2iQ(x))e 2i *. (80) 



We obtain 



24X 2 r^^Oftn {u k ( V i) 2 g* k (v) 2 ) = ^ 2 kM (81) 

where the amplitude M. is given by 



1 1 



M = -(Pcos2x + 2Qsin2x)Ci(-2x) + (Psin2:r-2Qcos2x)Si(-2x) -P+ — + - ) Q. 

X 2 \X d X J 

(82) 
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The cross terms in the Hamiltonian (1581) in Fourier space become 
lyflX f ,, / , A 1 \ 2v/6X /" d 3 A; 



d 3 x I 7r'o; H — na 
i] J \ i] 



V J (2vr) 3 



(83) 



where 



Mv) = u 'k(v) + -Mv)- 



(84) 



They do not affect (1691) and (17D|) at tree-level but are responsible for the remaining terms. 
After eliminating the contributions corresponding to disconnected diagrams, they yield 



H) + dZ2D + ([73D = -48X 2 (27r)^(k + p) 

/v rvi ]_ 
dm / drj 2 KMM^H^iK^i)^??) 2 + (c.c.)) (85) 

J vim 

+48X 2 (2vr) 3 5(k + P ) f ** [" ^v k ( Vl )u k ( Vl )v* k (rj 2 )ul( V2 )\g k ( V )\ 2 . 

J Vl J V2 

The integrals are not as bad as they look once the relations 



and 



v k (Vi)u* k (r]i) 



Vk(r]2)u k {rj2 



■"' f 1 ! % \ ,-,.!-,. !/,,,., 



2 a/cI V c s^7i 



"' / 1 _ ? 1 ^-i(l+c s )fc772 



2^ V C s^2 



(86) 



(87) 



are taken into account. The result is 



flUD + dZ2D + ([73D = 6X 2 £v4(27r) 3 5(k + p) 



(88) 



with the amplitude 



A 



4(3 + X) 



c 3 x 3 



Ci(— (1 + c 3 )x) sin(l + c s )x + Si(— (1 + c s )x) cos(l + c s )x 

Z, .1 



+^Q{x) 



Ci(— 2x) sin 2x + Si(— 2x) cos 2x 

2x 



— P{x) [Ci(-2a;) cos2x + Si(-2z) sin2x] 



■—z-Qix) [Ci( — (1 + Cg)x) cos(l + c s )x + Si(— (1 + c s )x) sin(l + c s )x\ 
cjx 

2 f x dy f y dz 
-— I — / — \P{x) sin(2x — y — z) — 2Q(x) cos(2x — y — z)\ sincJy — z) 
cjj y J z 



19) 



+— (P(x) sin2x - 2Q(x) cos2x) [Ci(-(1 - c s )x)Si(-(l + c s )x) - Si(— (1 - c,)s)Ci(-(l + c a )x)\ 



+— (P(x) cos2x + 2Q(x) sin2x) [Ci(— (1 - c s )x)Ci(-(l + c s )x) - Si(-(1 - c,)z)Si(-(l + c s )x)] 
cj 

If,, 2(3 + X 
+- [P(x) + -± =- 



2\2 



x- 



\ (Ci(-(l + c s )x) 2 + Si(-(l + c s )a;) 2 ) 
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The a-a correlation (164")) is needed only at the leading order and given by 

dMD = |j fc (77)| 2 (27r) 3 5(k + p). (90) 

Finally, the cross correlation flBol) is calculated as 

§B=mf dviHjim) (g k a k + ^o f _ k ) j* p bl p \0) + (h.c.) 
+i(0\ J" drnHriTu) (j k b k + f k bl^ ^ol p |0) + (h.c.) 

= W6Zi f —[vt(m)ut(Vi)9k(v)Jk(v)-(o.c.)] = 2y/6IkB (91) 

J m 



with B given by 



3 + X 2 / 3\ 3 / 3 + X : 



r2 



_ 1_ / 3(3 + X 2 ) r 3 + X 2 \ f _ _3_\ 1 

c s \ c s x 2 V x2 J V c s x2 / J 

x (Ci(— (1 + c 8 )x) cos(l + c s )x — Si(— (1 + c s )x) sin(l + c s )x) 



x (Si(— (1 + c s )x) cos(l + c s )x + Ci(— (1 + c 8 )x) sin 1 + c s )x) . 

Combining all the pieces, the curvature power spectrum is given by 

tt2 / 07-2 

\2|- \ _ a ^J l„./„\|2 , at2 U \A 1 °- L l„'_^M2 , «-t-2jl /I , 10 t2i 



(in|a(77) 2 |in> = ^-f r/ 4 Ug k ( V )\ 2 + 4X 2 A;.M + — ^MV + 6X 2 A^ + 12X 2 A;£) , (93) 

where we used a = —1/Hrj. In the standard calculation for single-field models, one evaluates 
this quantity at the horizon crossing since ( is conserved, even at non-linear order. Here 
in contrast, ( still evolves beyond the horizon scale and therefore we are going to evaluate 
its value at the end of inflation rj — rjf. The scales of cosmological interest should be well 
outside the horizon at the end of inflation, which means — krjf <C 1. In the corresponding 
limit x — y 0, the power spectrum reduces to 

<in|C*fo/)» -► j|^ (l + 18VQ1 2 (In |%|) 2 ) . (94) 

The disadvantage in this strategy is to introduce extra errors by neglecting time-variation 
of H and the slow-roll parameters over a number of Hubble times. We assume that these 
corrections are subdominant since they are suppressed by e# and the other small parameters. 
In fact, the leading correction coming from varying H is proportional to e# In \krj\, which can 
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be safely discarded compared to X 2 (ln \ki]\) 2 . It is difficult to estimate the effect of varying 
I as it is already a correction from the non-adiabatic evolution of £. It is expected to be 
proportional to 

(95) 



I __ e H Vh - V v 



HI 2 e H -e v 

multiplied by several powers of In \kr)\. Thus the approximation might break down when the 
e-folding number is too large. 



V. PRIMORDIAL GRAVITATIONAL WAVES AND PHENOMENOLOGICAL 
IMPLICATIONS 

A. Tensor Power Spectrum 



Let us now turn to the tensor perturbations 



7i j n "-ij , 



X% j - U i'j 



(96) 



where both h^ and Wy are trace- and diergence-free. Our quadratic action is 



f 



C t = -o- {KjKj - hij,khij,k) + jj^hijhij + — (uJi^ - uJij, k uJij,k) - c hiju'ij. (97) 



Normalizing the variables by 



it becomes 



C 



(2) 
T 



4/ 



Hij o/ijj, ojij jcjjj, 



o" 2c2 i i , 



(98) 



h' h' — h h 4- 

ij ij ll, tj,k'hj,k T| i „ 2 2 y -nj-nj 



1 



f" 



+ 2 I ^'A'j ~~ &ij,kOJij,k + -jCjijCjiA - — lu-j - — Wy J hij. 



2c 



r 



(99) 
(100) 



They are decomposed according to their polarization as 



hij{ri,x) = Y^ ^ijh'in,*), 4j4j = &s 



(101) 



s=l,2 



and a similar decomposition applies to Uij as well. Applying Fourier transform as before 

d 3 k 



h s (v,*) 



(2tt)* 



hi( V y^ + h^( V )e 



-ikx 



(102) 
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what one wants to compute here is the amplitude of the gravitational waves 

(27r) 3 <5(k + p)<m|/*fa) » = ^ E / (03 H (feW + ^kW) (%fa) + Wpfa)) l in >- 

(103) 
For this purpose, we only have to keep the leading order terms in Cj^ and therefore 

j ~ -2H, y ~ 2H 2 . (104) 

After discarding higher order corrections, the problem becomes a quantum field theory for 
the Lagrangian 



c = \ E [(#')')' ^ (v ^ )2 + ^ £ (^) a + ^)') 2 - < v **) a + ^( <D ') 2 1 ( 105 ) 

s=i,2 L ^ V . 

^-V^^H ((*■)' -V) &•. 

' 8=1,2 \ ' y 

One notices that the two polarization decouple from each other and the Lagrangian is of 
the same form as the scalar modes with different coupling constants and twice many fields. 
The interaction Hamiltonian for the each polarization mode is given by 

H} = -^ j d 3 x(h s ) 2 + 2V€H - ^ j d 3 x ((h s )'u} s + -kuA . (106) 

Hence, if we focus on the corrections up to first order in e# and e v , we only have to repeat 
the calculations in the previous section with replacements 



Cgp / (X fh 



a^/^^^-k + ^^lXk + ^XVl) for m (107) 

and 

2y/ * H ~ ^ J-0J3 (vkUkOtft* + v k u* k b£ai + v* k u k aib s k + ^afo-k) fOT ©, (108) 

and multiply the result by a factor of two to add up the polarizations. We also substitute 
itfe (77) for gk{ff)i which simplifies the calculations. In the end, one obtains 

(in|/i fe (77) 2 |m) = ^L (l + 1 - ^M + 2{e H - e v )M (109) 

with the amplitudes given by 

M = — ^- (Ci(-2x) cos 2x - Si(-2z) sin 2x) (1 10) 



X 2 X 2 
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— (Ci(-2x) sin 2x + Si(-2a;) cos 2x) , 
A= \ ~ x > S i n 2x + - cos 2x + —^- (Ci(-2a;) 2 + Si(-2x) 2 ) 



'MD 



+ 



l-x 2 



x- 



x/- 



In \x\ + 4 ) (Ci(-2x) cos 2a; + Si(-2x) sin 2x) 



+ — (In \x\ — 6) (Ci(— 2a;) sin 2x — Si(— 2a; ) cos 2a;) 
a; 

— cos(2a; — y — z) — sin(2a; — y — z) 



J y J z ix 

In the limit x — > 0, it becomes 



x- 



sin(y — z). 



AH 



m\h k (vf) 2 H ->• -TT i l + 4 ( e " - e ^) ( ln l%D ] • 



A- 3 



;ii2) 



B. Vector Mode 

Before discussing the observational implications, we shall briefly look at the vector per- 
turbation 



a% — v i-, X l j — K (i,j) + e ijk^k 



(113) 
(114) 



for completeness. As before, we take the flat slicing Fi = 0. The Lagrangian is given by 



f2 x2 f2 

£(2) _ ,J I I J , f2\/ \/ i \ \ 



^ijk^k,l^j,il H ~T^'i,j^ > i,j ® ^i I -^ ^i,jj ~i -^ ^ijk^j,k 



4 
2 



a 
1 



a 



a 



Varying S^ yields 



of 2 
V 2 ^ = — —A' (V 2 ^ + 2curlAi) 
a ' 



Similarly, Vi is non-dynamical and solved as 

1 



V 2 Vi = --V 2 /t- + curlA-. 



(115) 



116) 



117) 



Perhaps there are two promising choices of the gauge for the vector fields; A« = and Ki — 0. 
The action becomes 

£ = X (V«J) • (V<) - £ (V 2 k,) (V 2 ^) (118) 
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for the former and 

C = fKK - y (VAO ■ (VAO (119) 

for the latter. In either way, we have a free massless vector field with the propagation speed 

C S = 1/V2. 



C. Phenomenological Consequences 

If we assume instantaneous reheating whereby all the energy of the inflaton and gauge 
fields is dumped into a single relativistic fluid, the scalar curvature ( and the gravitational 
wave hij are conserved until re-entering inside the Hubble horizon and are observable through 
the CMB. The vector perturbation generated by the quantum fluctuation will quickly decay 
away and not be observed. For the scalar mode, the relevant parameter is the spectral tilt 
ns — 1. From ( 19^1) . this model predicts 

ns ~ 1 = ~±- In (^<C, 2 » = 36V f ln lkVfl 2 , (120) 

dink V ^ kn 1 + 18^X2(1^% |) 2 ' V ; 

which is negative for \kf]f\ <C 1, whence the spectrum is red. There is no contribution from 
the time- variation of H, €h, v and X since the spectrum was evaluated at the end of inflation, 
not the time each mode crossed the horizon. It can be seen that 

|n s -l| <min(-36v / 6X 2 ln|fc77 / |,- ^ J , (121) 

which means the spectral tilt doesn't necessarily impose a stringent constraint on the value of 
X. We note that ln \krjf\ is the e- folding number counted from horizon exit of the mode with 
wavenumber k until the end of inflation, whose value is model-dependent. If, for instance, 
we take ln \kr]f\ ~ —50, it yields 

|ns-l|~i (122) 

which is nicely consistent with WMAP 7-year [36J], even with X ~ 1. The running of ns can 
also be computed as 

5£*< n *- 1 > = £ji^-<'"- 1 >'' (123) 

which is safely small as long as ns — 1 is small. Similarly, the tilt of the tensor spectrum is 



given as 

7Zt 



J- In (e { ht)) = 8(^-e,)ln|%| (124 ) 

dink V kl> i + 4 (e H -e v )(ln|A;r // |) 2 
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Thus, the spectrum of gravitational wave is also red. However, it is much closer to scale 
invariance than that of scalar if X <C 1 . 

As can be seen from (j!12p . the amplitude of the tensor mode itself no longer provides 
the unambiguous information about the energy scale of inflation as it receives a potentially 
significant correction. The tensor-to-scalar ratio r can be used to determine e# and e v . 
Recalling X 2 = (e# — e ¥ ,)/e ¥ ,, it is given by 

m , i + 4(^)0^1)' 

(CI) H t v + isVe{(„-( v )(\n\k % \) 2 

Hence, the tensor to scalar ratio is suppressed compared to the single-field slow roll inflation. 
The suppression is stronger when the e-folding number is greater and the scalar kinetic energy 
is subdominant. 

In summary, we have seen that there are several different regimes that are consistent with 
the observations made so far, as far as the above three quantities are concerned. They are 
classified in the following. 

X 2 <C 1: Recalling the formula ( 123]) . this occurs when the energy density of the gauge 
fields is much smaller than the scalar kinetic energy density. If the e-folding number 
experienced by the modes at CMB scales is of order hundred or so, the model predicts 
a slightly red scalar power spectrum and almost scale invariant gravitational waves. 
Since it means e# ~ e^, the tensor-to-scalar ratio is unchanged from the ordinary 
single-field slow-roll inflation. 

X ~ 1: Although it means the dominant contribution to the scalar power spectrum comes 
from the terms proportional to X 2 , whose origin is the interaction between the inflaton 
and the gauge fields, the model still predicts an observationally consistent spectral 
tilt. On the other hand, even though the background energy density of the gauge 
fields is comparable to the scalar kinetic energy, the spectrum of gravitational waves 
is not very much affected for — In \kr)f\ ~ 50. The tensor-to-scalar ratio is suppressed. 
Of course, it should be noted that the perturbative approach cannot be trusted and 
non-linear effects may significantly modify the results. 

— In \krjf\ S> 50: Surprisingly, this regime is viable regardless of the value of X. The spectral 
tilt of scalar curvature is ~ — 1/ In \ki]f\ and the tensor mode is suppressed by a factor 
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of e 2 H with respect to the scalar mode, which makes it practically undetectable. A 
caveat is that time- variation of the background quantities (e.g. H and X), which was 
neglected in this article, might be important in this regime. 

VI. CONCLUSION 

In the present article, we investigated the linear perturbation of the inflationary model 
with a triad of background gauge fields coupled to the inflaton. We characterized the accel- 
erated expansion by introducing four parameters e#, e^, r/ H and r/^ which are generalizations 
of the usual slow-roll parameters. The general second order Lagrangian was derived and 
irreducible mode decomposition was carried out according to the transformation property 
under the spatial rotations, and the internal "rotation" of the triad of gauge fields. The 
scalar- and tensor-mode power spectra were computed by employing the in-in formalism on 
the de-Sitter background. We found the scalar fluctuation is characterized by the parameter 
X, which is potentially of order unity, while its tensor counterparts are smaller by an order of 
e v . The enhancement in the scalar correction is due to the steep gauge-kinetic function. The 
observational implications were studied by looking at the spectral tilt and tensor-to-scalar 
ratio. The generic prediction is that the spectra are red with the stronger effect for the 
scalar mode and the tensor-to-scalar ratio tends to be suppressed. The magnitudes of tilt 
and the suppression depend on X and the e-folding number. The structure of the corrections 
is such that even accurate measurements of those quantities are unable to impose strong 
constraints on X or the other small parameters because of the involvement of the e-folding 
number. As it stands, both weak and strong background gauge fields in the unit of inflaton 
kinetic energy are consistent with the results from WMAP. 

There are several implications for the preheating of auxiliary particles by the inflaton. 
The coupling needed to produce those particles back-reacts onto the scalar perturbation and 
can modify the power spectrum significantly. While gravitons can also be produced even 
in this isotropic background, which could affect the determination of the energy scale of 
inflation by measuring the amplitude of tensor mode, the effect is much smaller than that 
on the scalar mode, at least in the present model. 

Since it appears that the model passes the observational tests at linear order, it will be 



worth looking at the higher order corrections, namely non-Gaussianity. A recent work 29] 
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suggests that the gauge-kinetic coupling can lead to a strong signal, although the analysis 
was done for a vanishing background gauge field. In the above calculations, we have ignored 
the entropy modes and vector mode, which should affect the result at non-linear order. It 
is certainly interesting to carry out a thorough analysis taking into account all the modes 
involved, which should be possible for this isotropic model, and will be presented in the near 
future. 
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